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ABSTRACT 

We have obtained the best fit to the light profiles of a luminosity limited sample of 
elliptical and SO galaxies with a power law r 1 /™, letting the exponent remain free rather 
than keeping it fixed at 1/n = 1/4 as in the well known de Vaucouleurs formula. 
The introduction of a free parameter in the fitting formula (ranging from n = 0.5 for 
(r e ) = 0.3 kpc to n = 16 for (r e ) = 25 kpc) is justified by the existence of a good 
correlation between n and the global galaxian parameters, such as total luminosity and 
scale-radius. This result seems to be in line with the segregation of properties between 
the 'ordinary' and 'bright' families of early-type galaxies, and has consequence for the 
claimed independence of the shape of galaxy profiles with respect to the Fundamental 
Plane parameters. 
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1 INTRODUCTION 

Several attempts have been made in the past to find 'good' 
empirical representations of the light profiles of galax- 
ies, to be used as tools for quantitative classification and 
parametrization, for the identification of components and/or 
signatures of interactions, and as tests for models and predic- 
tions of numerical simulations (cf. the review by Capaccioli 
1989). 

Among the formulae proposed so far, the r 1 ^ law, in- 
troduced by de Vaucouleurs (1948) to represent the light 
profiles of ellipticals and of the spheroidal components of 
disc galaxies, has justifiably gained a popularity which has 
often been confused with 'universality'. 

The r 1 ^ formula 



I(r) 



'eexp{-*[(£) 1/4 -l]} 



(1) 



has 'scale' parameters only: a characteristic radius r e , and 
the corresponding surface brightness 7 e . The constant k 
is usually chosen in such a way that r e is the radius of 
the isophote containing half of the total luminosity; then 
k = 7.6692, and r e and 7 e are named effective radius and 
effective surface brightness. Useful formulae and tables re- 
lated to the r 1 ^ law may be found in Poveda, Iturriaga fc 
Orozoco (1960), de Vaucouleurs (1962), Young (1976), and 
Mellier fc Mathez (1987). 

The r 1 ^ formula has succeeded in reproducing, with 
a remarkable accuracy, the luminosity profiles of quite a 
few E galaxies. For instance, according to de Vaucouleurs 
fc Capaccioli (1979) and Capaccioli et al. (1990a), the r 1 ^ 



fit of the surface brightness distribution of the standard 
El NGC 3379 gives observed— calculated (O — C) residuals 
smaller than 0.08 mag over a 10-mag range. Also, the r 1 ^ 
law has been often used as a reference for the class of E 
galaxies in testing and interpreting the results of numeri- 
cal models of violent relaxation (van Albada 1982), merging 
of galaxies (Barnes 1988), and tidal stripping of ellipticals 
(Aguilar fc White 1986). While these models produce r 1 ^- 
like density profiles, a physical basis for this law is still lack- 
ing. 

Moreover, the analysis of the systematic deviations from 
the r 1 ^ behaviour has proven an effective tool for the dis- 
covery of secondary components (such as stellar discs), sub- 
structures (lenses, shells, ripples, etc.) and tidal effects seen 
as a shallowing or a steepening of the light profiles in the 
outer parts (cf. Kormendy 1977; Prugniel, Nieto fc Simien 
1987). Even the very first identification of the two fundamen- 
tal photometric components of disc galaxies (de Vaucouleurs 
1958, 1959) originated from the intuition that bulges of spi- 
rals might have r 1 ^ light profiles. 

Modern photometry of galaxies has shown, however, 
that for a large fraction of elliptical galaxies and of spheroids, 
the r 1 ^ law is only a first-order approximation: often it fits 
well only a limited interval of the luminosity profile. 

For instance, Michard (1985) noted that systematic de- 
viations from the r 1 ^ fit are similar for galaxies of similar 
luminosities. The best match seems to occur for objects with 
Mb = —21, with typical deviations of the order of 0.1-0.2 
mag arcsec -2 (Capaccioli 1985). 

Schombert (1986) found that, for a large fraction of the 
brightest cluster galaxies, the r 1 ^ interpolation is applicable 
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to the middle part of the light profile only (between 21 and 
25 B mag arcsec -2 ). More importantly, he pointed out that 
the shape of the light profiles is a unique function of the total 
luminosity: they become shallower in a (r 1 ^, fi) diagram as 
the luminosity increases (see his fig. 14). 

Binggeli fc Cameron (1991) have shown that, at the 
faint end of the galaxy luminosity function, dwarf ellipti- 
cals are better described by exponential or by King model 
profiles. By examining a wider range of luminosities, the 
same authors conclude that 'profile type is strongly corre- 
lated with the total magnitude of the galaxy'. 

As for bulges of SOs, it has been noted (e.g. Capaccioli 
1987, 1989) that the minor axis light profiles of edge-on ob- 
jects deviate systematically from the r 1 ^ trend [an upward 
concavity in the {r l ^ l ,fi) representation]. 

All of the above evidence casts some doubts on the uni- 
versality of the r 1 ^ law, and in particular on its reliabil- 
ity for deriving accurate effective parameters for early-type 
galaxies, as the results are too dependent on the portion 
of the light profile that is fitted^. This fact, and the uncer- 
tainty of the sky-background subtraction, is the main reason 
for the large discrepancies in the values of the effective pa- 
rameters of early-type galaxies reported by different authors 
(cf. Capaccioli, Caon fc D'Onofrio 1992b). 

Prompted by these considerations, we decided to study 
the behaviour of a generalized de Vaucouleurs formula (Ser- 
sic 1968), i.e. of a power law in which the exponent is a free 
parameter 1/n. Let us say explicitly that the emphasis of 
this paper is not in showing that the r 1 ^" law provides bet- 
ter fits than the simple de Vaucouleurs formula; improved 
representations of light profiles are just expected as an obvi- 
ous consequence of the introduction of a free parameter. Our 
aim is instead to investigate whether the parameter n has 
any physical justification by searching for correlations with 
global photometric parameters (effective radius and total 
luminosity). 

In Section §2 we give a brief account of the properties 
of the r 1/n law. A more complete description (luminosity 
density, dynamical characteristics and phase-space distribu- 
tion function) is given by Ciotti (1991). In Section §3 we 
present the application of the r 1 ^" law to the light profiles of 
a volume-limited sample of early-type galaxies in the Virgo 
cluster. Finally, in Section §4 we discuss the implication of 
the correlation between n and r e , also in relation to the 
dichotomy between the 'ordinary' and the 'bright' galaxy 
families claimed by Capaccioli, Caon fc D'Onofrio (1993). 



2 PHOTOMETRIC PROPERTIES OF THE r 1/n 
LAW 

In its simplest form the r 1 ^" law, originally proposed by 
Sersic (1968), is given by the formula 



7(7-) = 7 dex(-7- 1/ "). 



(2) 



^ We remark here that in principle r e is a model-independent 
parameter, being simply the radius of the isophote encircling half 
of the total luminosity. Therefore its identification with the scale- 
radius of the r 1 / 4 law is justified only for those galaxies that are 
perfectly fitted by the r 1 / 4 law over the whole light profile. 



Introducing the scale-radius re, it can be rewritten as 
J(t-) = J e dex j-6„ (I.y ,n -i J 
and, in logarithmic form, as 



l/r, 



1 



(3) 



(4) 



with c„ = 2.5 b n . The coefficient b n can be chosen in such 
a way that the scale-length r e is the radius encircling half 
of the total luminosity Li . The latter can be expressed as a 
function of r e and 7 e 



K n I e r. . 



(5) 



We have computed the coefficients b n and K n by numerical 
integration of an r 1 ^" galaxy model. In the range explored, 
0.5 < n < 16.5, they are very well approximated by the 
relations 



b n = 0.868 71 - 0.142 
and 

log(7f„) = 0.030[log(7i)] 2 + 0.441 log(Ti) + 1.079, 

both with a rms scatter of 0.014. 

For fitting purposes we will use the relation 

H(r) = A + Br 1/n 



(6) 
(7) 

(8) 



where, according to equation 4, the effective parameters of 
the j- 1 /™ model are 



-\- c n , 
(c„/.B) B . 



(9) 



We have added a superscript * to the scale parameters of 
the j- 1 /™ formula to prevent confusion with their model- 
independent analogues. 



3 AN APPLICATION OF THE 



LAW 



The r 1 /™ law has been fitted to the major and minor axis 
light-profiles and to the equivalent profiles of 52 early-type 
galaxies belonging to a sample of Virgo cluster members that 
is complete in luminosity: all E and non-barred SO galaxies 
of the cluster brighter than Bi = 14, or Mb = —17.3. The 
7?-band photometry that we use here is the result of the cou- 
pling of CCD images with deep Schmidt plates. This tech- 
nique (Capaccioli fc Caon 1989) provides accurate light pro- 
files spanning wide surface brightness ranges. More details 
about the properties of the sample and the data reduction 
can be found in Caon, Capaccioli fc Rampazzo (1990) and 
in Capaccioli et al. (1992a). 

The radial range of the light profiles over which the best 
fit with relation (8) has been computed goes from outside the 
region dominated by the seeing (typically a few arcsec, and 
fi b > 19.6 ± 1.1) out to where the uncertainty in the sky 
subtraction and some possible contamination from nearby 
objects render the light profiles unreliable. In other words, 
we fit the formula to the entire range of surface brightness 
covered by state-of-art photometric observations, removing 
only what is not measurable or just not directly measured. 
More specifically, deconvolution experiments (see below and 
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Caon et al. 1990) have been used to determine the boundary 
of the region where the surface brightness is changed by 
more than 5 per cent because of seeing convolution. As for 
the outer boundary of our fitting range, we have adopted the 
standard value of fiB = 27 mag arcsec -2 , where the average 
photometric error is not larger than 0.25 mag according to 
fig. 3 in Caon et al. (1990). This limit has been moderately 
varied to account for the properties intrinsic to each galaxian 
image. 

For some objects with peculiar inner light-profiles, the 
fit has been restricted to a smaller radial range. In a few 
cases, especially for very flattened discy galaxies, the fit of 
the r 1/n law was simply impossible. Seven galaxies out of 52 
could not be fitted at all because of the presence of promi- 
nent disc components, dust lanes or tidal distortions or sim- 
ply because of the irregular shape of the light profiles (see 
Table 1). For the remaining 45 galaxies, 12 major axis light 
profiles and five equivalent profiles could not be fitted by 
the r 1/n law. In some cases the fit required an additional 
(exponential?) component, while in others the profiles are 
intrinsically distorted. 

In order to minimize the possible repercussions of seeing 
blurring for results, we have repeated the fitting procedure 
with light profiles deconvolved according to the technique of 
Bendinelli (1991), and using in all cases a single Gaussian 
point spread function (PSF) with the stellar full width at 
half-maximum (FWHM) measured by Caon et al. (1990). 

In practice, for each galaxy profile we have scanned a 
fine grid of values of n, and computed the rms scatter of the 
(O — C) residuals. We have then chosen as the best value for 
the exponent the one giving the smallest scatter <Tbf. Note 
that for many objects the residuals are due to noise and 
do not reflect systematic departures from the r 1 ^" law. The 
resulting values of n range from 1 to 15. Two examples of 
the quality of the fits are shown in fig. 1, for a high and a 
low value of n. 

In order to estimate the uncertainty associated with the 
determination of n, for each galaxy (and for the major axis 
profile only) we computed the values n\ < n and 712 > n at 
which the rms scatter a of the (O — C) residuals exceeds by 
25 per cent the value of the best fit. For n ^ 10 we obtain 
typically 



An 



n 2 



n ~ 0.25ti. 



(10) 



For n > 10, the above relation still holds for n\ , while n^—n 
can reach values as large as n itself. For example, if a light 
profile is best fitted by an r 1 ^ law (typical scatter <Tbf ~ 
0.05), the use of an r 1 ^ 3 or an r 1 ^ 5 laws would cause a to 
increase by 25 per cent. On the other hand, when n > 8 is 
the best-fitting result, the same increment of a is achieved 
by using an r 1 ^ 8 ^ 2 ) law. 

For each galaxy of our sample, Table 2 provides the 
identification (column 1), the morphological type accord- 
ing to Binggeli, Sandage fc Tammann (1985) (column 2), 
the total luminosity in the B-band, the logarithm of the 
(model-independent) effective radius r e (in kpc), and the 
corresponding surface brightness 7 e (columns 3 to 5) from 
Caon et al. (1990) and Trevisani (1991), the borders //start 
and // e nd of the surface brightness range of the fit (columns 
6 and 7), the values of n for the major and minor axes (n ma j 
and Timm) and for the equivalent light-profile (?ieq) (column 
8), the rms scatter <Tbf of the fit (column 9), the coefficients 
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Figure 1. Observed major axis £?-band light profiles of the SO 
galaxy NGC 4552 and of the E3 NGC 4464 (filled circles), fitted 
by an r 1//n law (solidline) with n = 13.9 and n = 2.4 respectively. 
Crosses mark the borders of the fitting interval. The relative un- 
certainty of the data is indicated by the error bars, barely visible 
everywhere but in the outer parts. The dashed lines reproduce 
the r 1 / 4 laws best fitting the data within the same interval. The 
dotted and dot-dashed lines represent the results of the convolu- 
tion of the r 1//n and the r 1 / 4 models with a single Gaussian PSF 
and the proper FWHM (2 arcsec for both galaxies, according to 
Caon et al. 1990). 



